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ABSTRACT 

We use the projected correlation function w p (r p ) of a volume-limited subsample of the Sloan Digital Sky Sur- 
vey (SDSS) main galaxy redshift catalogue to measure the halo occupation distribution (HOD) of the galaxies 
of the sample. Simultaneously, we allow the cosmology to vary within cosmological constraints imposed by 
cosmic microwave background experiments in a ACDM model. We find that combining w p (r p ) for this sample 
alone with the observations by WMAP, ACBAR, CBI and VSA can provide one of the most precise techniques 
available to measure cosmological parameters. For a minimal flat six-parameter ACDM model with an HOD 
with three free parameters, we find Q, m = 0.278^;^, ct 8 = 0.812+g;^f, and H = 69.8^ kms" 1 Mpc" 1 ; these 
errors are significantly smaller than from CMB alone and similar to those obtained by combining CMB with 
the large-scale galaxy power spectrum assuming scale-independent bias. The corresponding HOD parame- 
ters describing the minimum halo mass and the normalization and cut-off of the satellite mean occupation are 
M min = (3.03!°;^) x lO 12 /?" 1 M Q , M\ = (4.58^;gj) x I0 n h~ l M Q , and K = 4.44^. These HOD parameters thus 
have small fractional uncertainty when cosmological parameters are allowed to vary within the range permit- 
ted by the data. When more parameters are added to the HOD model, the error bars on the HOD parameters 
increase because of degeneracies, but the error bars on the cosmological parameters do not increase greatly. 
Similar modeling for other galaxy samples could reduce the statistical errors on these results, while more thor- 
ough investigations of the cosmology dependence of nonlinear halo bias and halo mass functions are needed to 
eliminate remaining systematic uncertainties, which may be comparable to statistical uncertainties. 

Subject headings: cosmology: observations — cosmology: theory — galaxies: formation — galaxies: halos 



1. INTRODUCTION 

Over the last several years, halo occupation models of 
galaxy bias have led to substantial progress in characteriz- 
ing the relation between the distributions of galaxies and dark 
matter. Gravitational clustering of the dark matter determines 
the population of virialized dark matter halos, with essen- 
tially no dependence on the more complex physics of the sub- 
dominant baryon component. Galaxy formation physics de- 
termines the halo occupation distribution (HOD), which spec- 
ifies the probability P(N\M) that a halo of virial mass M con- 
tains N galaxies of a given type, together with any s patial 
and velocity biases of galaxies within halos ( Kauffmann et al. 
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; Benson et al. 2000; Berlind et al. 2003; Kravtsov et al. 


2004 


). Given cosmological parameters and a specified HOD, 



one can calculate any galaxy clustering statistic, on any scale, 
either by populating the halos of N-body s imulations (e.g., 
Uing et aljfl998l l2002t iBerlind & Weinberdl200l or by us- 
ing an increasin gly powerful array of analytic approximations 
(e.g..|Peacock & SmitlJ2000tlSeliakl2000tlScoccimarro et all 
l2001[|Takada & Jairl2003T see lCoorav & Shethl2002l forare- 
cent review). 

The 2dF Galaxy Redshift Survey (2dFGRS) dColless et al" 
1 20031) and the Sloan Digital Sky Su rvey (SDSS) JYork et al 
hood lAbazaiian etai]|2003L 12004b!) allow galaxy clustering 
measurements of unprecedented precision and det ail, making 
them id eal data sets for this kind of modeling. Zeha vi et alJ 
(2004a) (hereafter Z04a) show that the projected correlation 
function w p {r p ) of luminous (Mo.i r < -21) SDSS galaxies ex- 
hibits a statistically significant departure from a power law, 
and that a 2-parameter HOD model applied to the prevail- 
ing ACDM (cold dark matter with a cosmological constant) 
cosmology accounts naturally for this departure, reproducing 
the observed w p (r p ). Here, Mo .i r is the absolute magnitude in 
the redshifted r band, with observed magnitudes K-corrected 
to rest frame magnitudes for the SDSS bands blueshifted 
by z = 0.1, the median redshift of the s urvey ( B lanton et alJ 
I2003al) . iMagliocchetti & Porcianil i2003l) have applied a sim- 
ilar type of analysis to w p (r p ) for a fixed cosmology in the 
2dFGRS. The halo model plus HOD has also been used to 
suc cessfully describe the cluster ing of Lyman-break galax- 
ies JPorciani & Giavaliscoll2002l) . high-reds hift red galaxie s 
JZheng|l2004l) . as well as 2dF quasars JPorciani et al J l200l . 
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Recently, it was shown that large-scale overdensities are not 
correlated with galaxy color or star formation history at a 
fixed small-scale overdensity, supporting the HOD ansatz that 
a galaxy's properties are related on ly to the host halo m ass 
and not the large-scale environment (Blanton et al. 2004a). 

In this paper, we go beyond the Z04a analysis by bring- 
ing in additional cosmological constraints from cosmic mi- 
crowave background (CMB) measurements and allowing the 
HOD and cosmological parameters to vary simulta neously. 
This investigation complements that of Zehavi et al. (2004b) 
(hereafter Z04b), who examine the luminosity and color de- 
pendence of galaxy HOD parameters for a fixed cosmology. 
It also complements analyses that combine CMB data with 
the large-scale powe r spectrum measurements from the 2dF- 
GRS or SDSS _ (e.g. iPercival et alJl200ll ISoereelet alll2003l 
Tegm arket alJl2004l) . Such analyses use linear perturbation 
theory to predict the dark matter power spectrum, and they 
assume that galaxy bias is scale-independent in the linear 
regime. It also complements HOD and cosmological pa- 
rameter determin ation approaches u sing galaxy-galaxy lens- 
ing in the SDSS (Seliaket al. 2004a) and their combination 
with Lyman-q fores t clustering in the SDSS quasar sample 
(Selia k et alj|20 04b). Our analysis draws on data that ex- 
tend into the highly non-linear regime, and in place of scale- 
independent bias it adopts a parameterized form of the HOD 
motivated by theoretical studies of galaxy formation. 

2. THEORY 

We explore spatially flat, "vanilla" cosmological models 
that have six parameters, (£l/,/i 2 ,fl c h 2 ,0 s ,ln(A),n,r), where 
flt> and il c are fractions of the critical density in baryons 
and cold dark matter; Q s is the angular acoustic peak scale 
of the CMB, a useful proxy for the Hubble parameter, Hq = 
lOO/i km s" 1 Mpc" 1 ; A and n are the amplitude and tilt of the 
primordial scalar fluctuations; t is the optical depth due to 
reionization. 

In the halo model of galaxy clustering, the two-point cor- 
relation function of galaxies is composed of two parts, the 
1-halo term and the 2-halo term, £(r) = 1 + £u,(r) + ^,2h(f), 
which represent contributions by galaxy pairs from same ha- 
los and different halos which dominate at small scales and 
large scales, respectively. Here, the correlation function is cal- 
culated at the effective redshift of our observed SDSS sample 
at z = 0.1, which is a nontrivial modification since errors on 
the amplitude of the power spectrum at small scales (i.e., erg) 
are found be comparable to the growth factor shift at z = 0.1 
(see below ). The calculation of the 1- halo term is straightfor- 



ward (e.g- lBerhnd & Weinber320 02): 



dn (N(N-l)) M 



2irr 2 nj J dM 2 



' -F'i^r- I dM. (I) 



2R vh (M) \2R 



where n g is the mean number density of galaxies calculated 
from the HOD and halo model, dn/dM is the halo mass func- 
tion, (N(N- 1))m/2 is the average number of galaxy pairs in 
a halo of mass M, and F(r /2R v \ r ) is the c umulative rad ial dis- 
tribution of galaxy pairs (Berlind & Weinbe rg! 120021 IZhengl 
2004). 

For the 2-halo term term, in order to reach the accuracy 
needed to model the SDSS data, we include the nonlinear evo- 
lution of matter clustering and the halo exclusion effect. This 



is done in Fourier space, 
1 

2^2 



6h(r): 



where 



Pes.(k) — P mm (k) 



l g Jo 



M m 



P^ik^—dk, 



dn 



(2) 



dM — (N(M))b h (M)y g (k,M) 



(3) 

The mean occupation of halos of mass M is (N(M)), y g {k,M) 
is the normalized Fourier transform of the galaxy distribution 
profile in a halo of mass M. We approximate halo exclusion 
effects in two-halo correlation separations of r by choosing 
the upper limit of the integral in Eq. (|3j such that M ma x is the 
mass o f a halo with virial radius r/2, as incorporated in lZheng l 
(2004), Z04a, and Z04b. The importance of the nonlinear 
matter power spectrum and halo exclusion in accurately mod- 
eli ng the two-halo galaxy cor relation functio n was also fou nd 
by Maglio cchetti & Porcianil ll2003) and lWang et alJ d2004). 

In order to accurately include the dependence of the halo 
modeling of galaxy clustering for a varying cosmology, we 
include cosmologically general (within ACDM) forms of the 
nonlinear matter power spectrum, halo bias, halo mass func- 
tion, and dark matter halo concentration. We use the nonlin- 
ear matter spectrum P^(k) of Smith et al.'s (2003) halof it 
code, modified to utilize a numerically calculated transfer 
function from the Co de for Aniso tropics in the Microwave 
Backgr ound ("CAMB. iLewis et al J |2000T) . based on CMB- 
FAST (Seliak & Zaldarriaga 1996). We use halo bias fac- 
tors bh(M) determined in the high-resolution simulations of 
Seliak & Warren (2004), along with its given cosmological 
dependence, which provides a better fit (lower \ 2 ) to our 
observational dat a than halo bias models based on the pea k 
backgroun d split dSheth & Tormenll 19991: ISheth et alJ bOOl). 
We use the Jenkins et al. (2001) spherical overdensity of 180 
[SO(180), Eq. (B3)] halo mass function, and include in its in- 
terpretation of the definition of hal o mass the variation in th e 
virial overdensity with cosmology (Brvan & Normal! 1998) 



18tt 2 + 82x-39x 2 
l+x ' 



(4) 



and its effect in relating the LTenkins et alJ J22Q1 ) mass 
funct ion to va rying cosmologies (see, e.g., IWhitd 12001: 
Hu & Kravtsov 2003). Here, x = fl m (z)- 1. The variation 
of the virial overdensity also changes the halo exclusion scale 
of R v i r (M max ) used in Eq. (j3ji. Since our luminous SDSS sub- 
sampl e populates ha los of M > 1O 12 M , the breakdown of the 
Uenkins et all ( 120011) mass-function fit at M < 10 10 M Q is not 
important. 

We assume that the average spatial distribution of satellite 
galaxies within a halo follows a Na varro-Frenk- W hite (NFW) 
density profile of the dark matter (Navarro et al. 1996), moti- 
vated by hydrodyna mic simulation results fWhite et all2001t 
IBerlind et alJ 120 03) and N-body simulation galaxy cluster- 
ing predictions with halos populate d by semi-analytic m odels 
(Kauffmann£^lJ1227yKauffman n~et all 199 9: Benso net alJ 
2000; Somerville et al. 2001 ). However, as a test, we drop this 
assumption of no spatial bias within halos between galaxies 
and dark matter and find it is not important (see |4] below). 
In the case of no spatial bias, each halo is assumed to have a 
cosmologically-dependent concentration 
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(5) 
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where 



1000 



0=11 



a2 m 
03 
0.05, 



1-1.7 



-1.6 



(6) 
(7) 



as found in fits to numerical resu lts for varying cosmologies 
bv lHuffenberger & Seliald ( I2003I Here. 



"eff = 



dinPUk) 



dink 



(8) 



where k* is the nonlinear scale such that Aj. n (fc*) = 1, M is the 
virial mass of the halo and M* is the nonlinear mass scale. 
There is a scatter about any mean concentration value and 
this could change the prediction of the shape of a given halo. 
However, as we describe below, our results are largely insen- 
sitive to the exact form of the concentration of the galaxies 
with respect to the dark matter. 

Our HOD parameterization for a luminosity-threshold 
galaxy sample (Mm, < -21 in this paper) is motivated by 
results of subst ructures from high-resolution dissipationless 
simulations of lKravtsov et alJ (120041) . The HOD has a simple 
form when separated into central and satellite galaxies. The 
mean occupation number of central galaxies is modeled as a 
step function at some minimum mass, smoothed by a comple- 
mentary error function such that 



(N cen (M)) = ^Erfc 



ln(M min /M) 



(9) 



to account for scatter in the relation b etween the a dopted 
magnitude limit and the halo mass limit JZhengetalJ l2004). 
(Note that the number of central galaxies is always N cen = 
or A^ C en = 1 by definition.) The occupation number of satellite 
galaxies is well approximated by a Poisson distribution with 
the mean following a power law, 



(N sat (M)) : 



M-nM„ 
Mi 



(10) 



where we introduce a smooth cut-off of the average satel- 
lite number at a multip le k > 1 of the m inimum halo mass. 
iGuzik & Seliakl ( 120021) and lBerlind et alJ d2003| ) using semi- 
analytic model calculations, and Kravts ov et alJ 112004 ). using 
high resolution A^-body simulations, found awl. The gen- 
eral HOD above is characterized by five quantities: M ra j n , M\, 
<7 cen , k, and a, and we refer to this as the 5p model. It pro- 
vides an excellent fit to predicti ons of semi-analyt ic models 
and hydrodynamic simulations (Zheng et al. 2004), in addi- 
tion to describing sub halo populations in N-body simulations 
JKravtsovet al 120041). 

3. OBSERVATIONS 

The SDSS uses a suit e of specialized instruments and 
data reduction pipelines JGunn et alJll998t iHogg et al]l2001t 
lPieretal.ll27)03t lStoughtonetalJl2002l) to image the sk v in 
five passbands (Fukugita et all 1 19961 [Smith et al .1 120021) and 
obtain spectra of well defined samples of galaxi e s and quasars 
(|Ei^enste2n^^lJj2p(^j3^ c h arc l s et alJ 120021 IStrauss et al.l 
120021 IBlanton et al. 2003b). For our analysis, we use Z04b's 
measurement of the projected correlation function w p (r p ) of 
a volume-limited sample of galaxies with Mo.\ r < —21. This 
sample is in turn selected from a well characterized subset 
of the main galaxy sample as of July, 2002, known as Large 
Scale Structure (LSS) samplel2, which includes -200,000 
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FIG. 1 . — Shown are the projected correlation function w p (r p ) of Mo.i r < 
—21 galaxies from the SDSS LSS samplel2 (points with lcr diagonal er- 
rors) and the best fit three parameter HOD model (solid). Only points with 
r < 20ft" 1 Mpc are used in the fit. Also shown are predicted w p (r p ) models 
with Q c (dotted) and erg (dashed) at ±3cr from their best fit values. HOD 
parameters and other cosmological parameters are held fixed. As seen here, 
the sensitivity to erg and Q c comes from both the amplitude and combined 
shape of the 1-halo and 2-halo regimes of w p (r p ). 



galaxies over 2500deg 2 of sky ( Bla nton et al. 2004b). We use 
the Mo.i r < -21 sample with the full w p (r p ) data covariance 
matrix from the jackknife estimates of Z04b. There are 26,015 
galaxies in the Mo.i r < -21 sample. 

The observed projected correlation function is obtained 
from the 2-d correlation function £(r p , 71") by integrating along 
the line of sight in redshift space: 



(11) 



w p (r p ) = 2 / £(r p ,ir)dir, 



where r p and it are separations transverse and parallel to the 
line of sight. We adopt 7r max = 40/r'Mpc (in the measurement 
and modeling), large enough to include nearly all correlated 
pairs and thus minimize redshift-space distortion while keep- 
ing background noise from uncorrected pairs low. Because 
our sample is volume-limited, we are measuring the cluster- 
ing of a homogeneous population of galaxies throughout the 
survey volume, which greatly simplifies HOD modeling. Fur- 
ther details of the sample and measurement are given in Z04b. 
In our analysis, we use 11 w p (r p ) data points in the range 
0.1/i _1 Mpc < r p < 20/i _1 Mpc, sufficiently below the projec- 
tion scale 7r max to avoid contamination of redshift space dis- 
tortions, though we have found that including points up to 
r w 40/i _1 Mpc, which have low statistical weight, does not 
alter our results. 

We also require our models to reproduce the measured 
mean comoving space density of our sample, n° bs = 1.17 x 

10~ 3 /! 3 Mpc~ 3 . This quantity has an uncertainty due to sample 

variance that can be written as cr s /n° bs = yj (<5|), where (<5 2 ) 

is the variance of the galaxy overdensity. We estimate (i5 2 ) by 
integrating the two-point correlation function over the volume 
of the SDSS Moa,- < -21 sample. To compute this integral 
we generate a large number of independent random pairs of 
points within the sample volume and sum £(r) over all these 
pairs. We use the Z04b correlation function and extend it to 
larger scales with the linear theory correlation function mul- 
tiplied by b 2 , where b = 1.4 is the large-scale bias factor for 
Mo.i r < -21 galaxies. We vary the number of random pairs 
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FIG. 2. — Shown are the marginalized posterior likelihoods for the cosmo- 
logical parameters determined from CMB+w p (r p )[Mo.i r < -21] with a three- 
parameter HOD in solid blue, that for CMB+SDSS 3D P g (k) in red dashed, 
and that for the CMB alone in yellow (gray) shaded. 



used and find that the integral converges at 10 5 pairs. Us- 
ing two million pairs, we find that the number density uncer- 
tainty due to sample variance is er fl /n° bs = 0.0377. There is 
also a shot noise Poisson uncertainty in the number density 
that, for this number of galaxies, is 0.0062. We add these 
two components in quadrature to obtain a total uncertainty of 
a? fflj. 



nf s = 0.0382, and therefore 



n° bs = (1.17 ±0.05) x 10" 3 /i 3 Mpc" 



4. RESULTS 



(12) 



For a given cosmology and HOD parameter choice, we use 
the predicted w p (r p ) to calculate the likelihood to observe the 
Mo.i r < -21 sample's w p (r p ) and n° bs . We combine this like- 
lihood with that for the model's prediction for the cosmic 
microwave background anisotropy temperature correlation 
and temperature-polarization cross-correlation to produce the 
WMAP (first year), ACBAR (£ > 80 0), CBI (600 < I < 
2000) and VSA (£ > 600) observations Jffinshaw et alfcOOl 
[ Verde et al J 120031 iKoeut et alJ 120031 iDickinson et alJ 120041 
Readhead et al. 2004»: We vary the six parameters for the 
"vanilla" ACDM cosmological model plus the five HOD pa- 
rameters: p = (£lt>h 2 ,£l c h 2 , @.s, ln(A), n, T,M m [ n ,M\ , k, a cen , a). 
The ranges allowed in our Markov Chain Monte Carlo 
(MCMC) sampling of parameters are chosen to avoid any ar- 
tificial cut-off of the likelihood space and are 

0.005 <n b h 2 < 0.1 
0.01 < fl c h 2 < 0.99 

o.oo5 < e s <o.i 

-0.68 < ln(A) < 0.62 
0.5 < n <1.5 
0.0K r <0.8 



10 9 M Q <M min < 5 x 1O 13 M 
10 10 M Q < Mi < 4 x 1O 14 M 
1 < k < 30 

< CTcen < 10 

0.2 < a <2.5. 



(13) 



Here, A is related to the amplitude of curvature fluctuations 
at horizon crossing, \AR\ 2 = 2.95 x 10~ 9 A at the scale k = 
0.05 Mpc. The angular acoustic peak scale Q s is the ratio of 
the sound horizon at last scattering to that of the angular diam- 
eter d istance to the surface of last scattering (Kosowskv et al. 
2002). 

To measure the likelihood space allowed by the data, we 
use a Metropolis MCMC method with a modified version 
of the Lewis & Bridld J2002I) CosmoMC code. We use the 
WMAP team's code to calculate the WMAP first-year ob- 
servations' likelihood, and CosmoMC to calculate that for 
ACBAR, CBI and VSA. After burn-in, the chains typically 
sample 10 5 points, and convergence and likelihood statistics 
are calculated from these. Since it is not known a priori which 
HOD parameters are most constrained by the w p (r p ) measure- 
ment, we use the Akaike and Bayesian Information Criteria 
(AIC and BIC) to determine w hich parameter s are statistically 
relevan t to describin g w p (r p ) JAkaikd[l974t ISchwarz|ll978l 
see also Liddle 2004). More parameters might well be needed 
once we have more data to constrain the HOD, but w p (r p ) 
alone doesn't provide enough information to demand it. 

Likelihood analyses were performed for several cases 
where some parameters were kept free and others were fixed 
to a physical limit, i.e. where the scatter in the mass- 
luminosity relation is unimportant (<7 cen = 0), the cut-off scale 
of the satellite galaxies is exactly that of the minimum mass 
M m j n (k = 1), or to a value (a = 1) predicted in the numerical 
simulations and semi-analytic mod els of satellite halo dis tri- 
butions in Berli nd et all (120031) and lKravtsov et alJ J2004I) . If 
we adopt all three of these constraints and allow only two pa- 
rameters, M m i„ and Mi, to vary to fit w p (r p ) and n° bs , then we 
obtain a poor fit. This model is an inadequate description of 
the data according to the information criteria (ABIC = 7.2 and 
AAIC = 12.5) relative to the three-parameter M m i n , Mi, and 
k model (3p). We also investigated a four parameter model 
(4p), varying M m i„, Mi, n, and cr cen with a = 1, as well as a 
five parameter model (5 p) varying all parameters in this HOD. 
Relative to the 3p model, the Ap and 5p models introduce 
new parameters that are not justified by the information crite- 
ria (ABIC > 6, cf. Tabled, since these models add freedom 
but yield only a very small reduction in x 2 - 

To assess the importance of one aspect of the halo model- 
ing, we performed a test on the 3p model whereby the NFW 
concentration cq [Eq. Q] of dark matter is replaced by that 
for the galaxies, c^', and is also left free in the MCMC within 

0.01 < Cg al < 200 and independent of the dark matter concen- 
tration of the halos. We find that the derived cosmological 
parameters and their uncertainties remain nearly unchanged 
from a model with no spatial bias, and the constraints on 
the galaxy concentration are consistent with no spatial bias: 
Cq" 1 ' = 11. I*!? 3. The marginalized values of the HOD parame- 
ters M min and Mj remain unchanged with varying Cq , though 
the error on the cut-off scale of the satellite galaxies k in- 
creases (k = 4.71+j 4j). This increase is expected since it is 
precisely the central distribution of satellite galaxies that is 
positively correlated to the one-halo galaxy distribution con- 
centration, with a correlation coefficient of r = 0.72. As a large 
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TABLE 1 

COSMOLOGICALPLUS HOD PARAMETERS, MARGINALIZED CONSTRAINTS WITH 68.3% C.L. ERRORS. 
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79 I 
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7 
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7tv J 
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n 1 t ^ n+0 0039 
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0.1107^00050 


0.1088^0058 


Oil 7A^' ul ' uy 
U - A1/D -0.0069 
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O.l04+O;0|3 




U - UZZ ~' -0.00080 


A rm/1-7+0 00084 
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FIG. 3. — Plotted are the Ax 2 < 1 range from the best fit for the HOD 
drawn from the MCMC chains for the 3p model in red (dark gray), 4p model 
in green (medium gray), and 5p model in yellow (light gray). 



Cy J makes the distribution of galaxies inside halos more con- 
centrated, to maintain the small-scale clustering, k increases 
to allow relatively more galaxies to be put in halos with larger 
virial radii and lower concentrations. 

Figure 1 illustrates the way that w p {r p ) constrains cosmo- 
logical parameters. Data points show the Z04b measure- 
ments, and the solid line shows the prediction of the best-fit 3p 
model. Dashed curves show the prediction of the w p (r p ) after 
as is perturbed by ±3<r relative to its best-fit value given in 
column 2 of Table 1, with all other cosmological parameters 
(and therefore the shape of the linear matter power spectrum) 
as well as the HOD parameters held fixed. Dotted curves show 
the prediction of w p (r p ) after changing O c , and thus the shape 
of the transfer function in the matter power spectrum, by ±3cr, 
with all other cosmological and HOD parameters fixed. The 
strength of the constraints derived from w p {r p ) stems from 
the combined relative dependence of the 1-halo and 2-halo 
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FIG. 4. — Shown are the marginalized 68.3% and 95.4% C.L. contours 
in ct 8 vs. U m for the WMAP+ACBAR+CBI+VSA (CMB) data alone (gray 
shaded), from the CMB + SDSS 3D P g (k) (orange/light-gray lines) and CMB 
+ SDSS w p (r p ) (blue/dark-gray lines) from the 3p HOD analysis presented 
here. 



regimes and therefore the overall shape of w p (r p ). 

The cosmological parameters' marginalized posterior like- 
lihoods for the 3p model are shown in Figure [2] Also shown 
for comparison are the marginalized likelihoods for the CMB 
plus SDSS 3D P g (k) [updated from Teg mark et aflfe004l) with 
new CMB results], and that from the CMB data alone. We 
also combine the CMB+w p (r p )(3p) measurement with the 
SDSS 3D P g (k) for a joint constraint on cosmological param- 
eters. Since the P g (k) data points included in the analysis are 
at wavelengths \ = 2-k jk> 30/r'Mpc, the information they 
contain is largely independent of that in the w p (r p ) data points 
at r p < 20h~ l Mpc. All parameters' best fit values and errors 
are listed in Tabled The resulting range of the HOD mea- 
sured for all models here are shown in Figure|3] 

Two-dimensional contours of f2 m and are shown in Fig- 
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FIG. 5. — Plotted are the 68.3% and 95.4% C.L. contours for the marginal- 
ized likelihoods for f!,„ vs. M m - m . The strong dege neracy (cor relation 
of r = 0.92) r oughly follows Q„, oc M mm , as expected I Zheng et aUI2002l 
IRozo et all2D 04) 



ure|4]and are compared to those obtained using CMB alone 
or CMB + P g (k). The anticorrelation of Q, m and a% from the 
w p (r p ) constraint seen in Fig.|4]arises from the anticorrelated 
degeneracy in these parameters in the one-halo component 
due to its dependence on the halo mass function which needs 
to maintain its amplitude at high halo masses, and the Omega,,, 
and a% anticorrelation in the two-halo component due to the 
amplitude-shape degeneracy of the dark matter power spec- 
trum (or dark matter correlation function). 

Important results to note from the Figures and Table are the 
following. Cosmological constraints obtained using CMB and 
Wp(fp) are substantially tighter than those from CMB alone, 
and they are similar in value and tightness to those obtained 
from CMB + P g (k) despite the introduction of new parame- 
ters to represent the HOD. The a% constraints using w p (r p ) 
are tighter than those using P P (fc); note that the latter estimate 
has dropped relative to that of Tegmark et al] ( 120041) because 
of the smaller scale CMB data. If we incorporate P g {k) con- 
straints in addition to w p (r p ), then parameter values change 
by less than I— a and error bars improve slightly. Our cos- 
mological parameter results also agree, within errors, with 
the recent results fro m SDSS galaxy bias and Lyman-a forest 
(Seliaket al. 2004b). The HOD parameters are partially de- 
generate among themselves, so adding parameters to the HOD 
model worsens the constraint on any one of them. However, 
within the range of models examined here, adding parameters 
to the HOD only slightly increases the error bars on cosmo- 
logical parameters. 

Since the small scales of the primordial power spectrum 
probed by w p {r p ) could be useful in constraining any devi- 
ations from a simple power-law primordial spectrum as well 
as a model including the suppression in power spectrum and 
mass function due to the presence of massive neutrinos, we 
performed an MCMC analysis including a running of the 
spectrum dn/dlnk about the scale k = 0.05 Mpc for the 3p 
HOD model as well a model including massive neutrinos. 
We find little evidence for runnin g, dn/dlnk = - 0.062^°^, 
comparable to the results of Sperg el et al.l ( 120031) . The halo 
model in the presence of massive neutrinos is applied such 



that the halo mass function, halo bias and halo profile is that of 
the cold dark matter alone since neutrino clusteri ng is a very 
small effect on these quantities ( Abazai ian et al.l2004al) . The 
presence of massive neutrinos is constrained to m v < 0.27 eV 
(95% C.L.) for each of 3 neutrinos with degenerate mass. The 
statistical errors from our analysis of the CMB plus this w p (r p ) 
measurement on dn/dlnk and m v are comparable to those 
from other cosmological parameter analyses, being smaller 
than those from the sh ape of the SDSS 3D P g (k) plus WMAP 
l lTegmark et al.ll2004l) . comparable to the WM AP plus 2dF- 
GRS 3D P g (k) plus modeled bias constraints ( Sperg el et alJ 
l200l . but not as stringent as those from modeling the galaxy 
bias in the SDSS from galaxy-gala xy lensing and cluste ring 
of the Lyman-a forest in the SDSS dSeliak et all2004albl) . 

5. DISCUSSION 

The remaining uncertainties in cosmological parameters in- 
troduce relatively little uncertainty in the HOD parameters, 
i.e., we now know the underlying cosmology with sufficient 
precision to pin down the relation between galaxies and mass. 
The strongest expected degeneracy is between the value of 
Sl m and the mass scale parameters M m ; n and M\, since one 
can compensate a uniform increase in halo mas ses by sim- 
ply shifting galaxies in to more massive halos llZheng et alJ 
120021 IRozo et al . 2004). The error contours for Sl„, vs. M m j n 
are shown in Figure [5] The degeneracy between these pa- 
rameters is strong, with a correlation of r = 0.96. While this 
degeneracy would cause large uncertainties in the values of 
M m i n and fi m if we used the galaxy clustering data alone, the 
combination of CMB and w p (r p ) data constrains Q. m fairly 
tightly, leaving limited room to vary the mass scale param- 
eters. Incorporating SDSS clustering measures that are di- 
rectly sens itive to halo m asses, such as redshift-space distor- 
tions llZehavi et al. 2002) an d galaxy-galaxy len sing measure- 
ments ( Sheldon et al. 2004; lSeliaketalJ l2004a). may further 
improve the il m con straints. 

As discussed by Berlind & Weinberg (2002), the galaxy 
correlation function places important constraints on HOD pa- 
rameters, but it still allows tradeoffs between different fea- 
tures of P(N\M) and (to a lesser degree) between P(N\M) 
and the assumed spatial bias of galaxies within halos. Ad- 
ditional clustering statistics such as the group multiplicity 
function, higher order correlation functions, and void prob- 
abilities impose complementary constraints that can break 
these degeneracies. Our analysis should thus be seen as a 
first step in a broader program of combining galaxy clus- 
tering measurements from the SDSS and other surveys with 
other cosmological observables to derive simultaneous con- 
straints on cosmolog i cal par ameters and t he gal a xy HOD [see 
Berlind & Weinberg] J200l : IWeinberd J2001: IZheng et al] 
(2002) for further discussion]. Van den Bosch et al. (2003) 
have been carrying out a similar program using the closely re- 
lated conditional luminosity function (CLF) method applied 
to the 2dFGRS luminosity and correlation functions (see also 
Ivan de n Bosch et al. 2004|- They find cr 8 = 0.78 ± 0.12 and 
O m = 0.25±g 07 in their analysis combined with CMB data, 
with both errors at 95% C.L. as given in that work. Our results 
are in agreement, within errors, with their determinations of 
erg and Q m . 

Besides the statistical error bars, there are two main sources 
of systematic uncertainty in our cosmological parameter esti- 
mates. The first is the possibility that our HOD parameteriza- 
tion does not have enough freedom to describe the real galaxy 
HOD, and that we are artificially shrinking the cosmologi- 
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cal error bars by adopting a restrictive theoretical prior in our 
galaxy bias model. For the 3p model, this is arguably the 
case, since it effectively assumes perfect correlation between 
the mass of a halo and the luminosity of its central galaxy. 
However, our 5p model is able to give an essentially perfect 
description of the predictions of semi-analyti c galaxy for ma- 
tion models and hydrodyn amic simulations (Zheng 2004; see 
alsolGuzik & Se liak 20021 iBerlind et alJ2003[ Kravtsov et al. 
120041) . so there is good reason to think that the error bars 
quoted for this case are conservative. This model still makes 
the assumption that satellite galaxies have no spatial bias with 
respect to dark matter within halos, but the concentration test 
in |4] shows that dropping this assumption has minimal im- 
pact on cosmological conclusions. In place of an HOD model, 
traditional analyses based on the large-scale galaxy power 
spectrum assume that the galaxy power spectrum is a scale- 
independent multiple of the linear matter power spectrum, so 
that their shapes are identical. Scale-independen ce in the lin- 
ear regime is expected on fairly general gr ounds (Cole sll993i 
[ Fry & Gaztanagal [T99l [W einberg fl995t iMann et alJ Il998l 
Ischerrer & Weinberg! Il998t iNaravanan et alJ 120001) . How- 
ever, it is not clear just how well this approximation holds over 
the full range of scales used in the power spectrum analyses, 
so although our HOD models are considerably more complex 
than linear bias models, our approach is arguably no more de- 
pendent on theoretical priors. In future work, we can use the 
HOD modeling to calculate any expected scale-dependence 
of the power spectrum bias, thus improving the accuracy of 
the power spectrum analyses and allowing them to extend to 
smaller scales. 

The second source of systematic uncertainty is the possibil- 
ity that our approximation for calculating w p (r p ) for a given 
cosmology and HOD is inaccurate in some regions of our pa- 
rameter space. The ingredients of this approximation have 
been calibrated or tested on N-body simulations of cosmolog- 
ical models similar to the best fitting models found here, so 
we do not expect large inaccuracies. However, there are sev- 
eral elements of the halo model calculation that could be in- 
accurate or cosmology dependent at the 10% leve l that is now 
of interest, including departures from the lJenkins et alJ d200 11) 
mass function, scale dependence of halo bias, and details of 
halo exclusion. Uncertainties in the halo mass-concentration 
relation and the impact of scatter in halo concentrations come 
in at a similar level, though the test in &\ again indicates 
that these uncertainties mainly affect the details of the derived 
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HOD, not the cosmological parameter determinations. With- 
out a comprehensive numerical study of these issues, it is dif- 
ficult to assess how large the systematic effects on our param- 
eter determinations could be, but we would not be surprised 
to find that they are comparable to our statistical errors. We 
plan to carry out such a study to remove this s ource of system- 
atic un cert ainty in future work; the papers of Selia k & Warrenl 
J2004I) and lTinker et alJ d2004) present steps along this path. 

Analyses of multiple classes of galaxies will allow consis- 
tency checks on any cosmological conclusions, since differ- 
ent classes will have different HODs but should yield consis- 
tent cosmological constraints. By drawing on the full range 
of galaxy clustering measurements, joint studies of galaxy 
bias and cosmological parameters will sharpen our tests of the 
leading theories of galaxy formation and the leading cosmo- 
logical model. With this current analysis alone, we find that 
the combination of CMB anisotropics and small-scale galaxy 
clustering measurements provides, simultaneously, tight con- 
straints on the occupation statstics of galaxies in dark matter 
halos, and some of the best available constraints on funda- 
mental cosmological parameters. 
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